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In this work assuming valid the equipartition theorem and using the normahzed q-expectation 
value, we obtain, until first order approximation, the hydrodynamics equation for the generalized 
statistics. This equations are different from those obtained in the context of the Boltzmann-Gibbs 
statistics. This difference is that now appears two transport coefficient that depend on the q-value. 
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I. INTRODUCTION 

The Generalized Statistics was first proposed by Tsal- 
lis ^ in 1988. It has attracted a lot of attention lately 
for it has been well applied to many systems, like fully 
developed turbulence ||] and self-graviting systems [||. 
For a periodically updated list please see reference Q . 
This generalization is based on the new entropic form, 
namely: 
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where fcs is the Boltzmann constant, pi is the proba- 
bility that the system is in a state labeled by i, g is a 
real parameter and the sum runs over all the allowed 
states. This entropy has the essential feature that it is 
non-extensive (non-additive). For two independent sys- 
tems, i.e., for two systems which have independent prob- 
abilities: = p(^)p(-^\ the entropy of the whole 
system is written as: 



(2) 



At that time it was used to calculate the physical quanti- 
ties average values without proper normalization. Later 
Tsallis, Mendes and Plastino Q have introduced a new 
form of calculating expectation values, the so called nor- 
malized q-expectation value. Let O be some physical ob- 
servable, its expectation value, in this context, is given 
by: 



< o>,= 



where: 



= [l-(l-<z)/3(i^,-^7,)]— 
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is the canonical ensemble distribution function, with Uq 
being the q-average internal energy known, a priori, from 
experiment. 

Tsallis introduced a cut-off condition ^ in this formal- 
ism. It says that if the argument of the distribution 



function is negative, this distribution vanishes, in other 
words, if 1 - (1 - q)P(Hi - Uq) < 0, = 0. 
In 1999, Boghosian ||] has obtained the hydrodynamic 
equations for the generalized statistics. He calculated 
them using the unnormalized q-expectation value and the 
well known Chapman-Enskog expansion. He showed that 
the equation of state and the internal energy for the ideal 
gas are q-dependent. He showed also that, in first order 
approximation, the conservation equations don't change 
and, in second order approximation, only the energy con- 
servation equation change. 

The first change is that the thermal conductivity coef- 
ficient is now q-dependent. It is given by, in the three- 
dimensional case: 



Aq + {1- q) 3mP 



where: 



Aq = l + il-q)- 
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T is the expansion parameter, p is the mass density of 
the system and f3 is the inverse of the temperature. 
There was, also, the appereance of a new term in this 
equation, a gradient of pressure, with a new transport 
coefficient, which was called anomalous transport coeffi- 
cient: 
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Obviously this coefficient vanishes in the extensive limit 
q^l. ^ 

We propose in this paper to obtain the hydrodynamic 
equations using a different method and using the nor- 
malized q-expectation value. We'll see that these new 
features will led us to the correct equation of state for 
the ideal gas and to different dependence on q of the 
transport coefficient. In section two we obtain the Con- 
servation Theorems for the conserved quantities in molec- 
ular collisions, we make the zero order approximation, 
that corresponds to the first order approximation in 
Boghosian Q calculations and we present the first order 
approximation (Boghosian second order approximation). 
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Finally we presents our results and conclusion in section 
three. 



II. NORMALIZED HYDRODYNAMIC 
EQUATIONS 

A. Conservation Theorems 

Here we present our way to obtain the hydrodynamic 
equations. We'll follow reference and assume the en- 
ergy equipartition theorem to be valid Q . We are study- 
ing an ideal gas in the 6-dimensional phase space. The 
one particle hamiltonian for this system reads: 
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where Ui — Vi — Ui, with —< Vi > the average velocity 
and < Ui >= 0, m stands for the mass of the molecule 
and the index i runs from 1 to 3 and indicates the com- 
ponents of the velocity. 

It consists in obtaining the equations which give the 
transport of the molecular quantities that are conserved 
in a molecular collision. To do this, we assume first that 
we are dealing only with binary collision and that the two 
colliding molecules are labeled 1 and 2. We are studying 
only the molecule 1. We also assume that the molecular 
chaos is still valid. 

To obtain the conservation theorems we consider the 
Boltzmann transport equation: 



d Fi 
— + v^di + —dy, 
ot m 



d\2d?Vf,d\fj\Vf - V,)\Tf,\^ iFf.,Ff, - F^F,) (8) 



where vi are the components of the thermal velocity of 
the molecules, Fi are the components of the external 
forces, which are considered here velocity independent, 
the index 2,/i,/2 refers to the second molecule before 
the collision and to the both molecules after the colli- 
sion, the 4-dimensional delta function and the transition 
matrix elements Tfi are related to the differential scat- 
tering cross section of the collision and the term between 
parenthesis is the molecular chaos assumption. Fi is the 
particle density which is the solution for (8) and is un- 
known. 

The equilibrium distribution function is given by the es- 
cort distribution: 



Fq = n- 
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where /3 is the inverse of the temperature /? — . The 
factor n is the particle density. It appears here because 



this distribution is described in the one particle phase 
space, so it is not simply a probability distribution, it is 
the number of particles in the system which have posi- 
tion between x and x + dx and velocity between v and 
V + dv. This function is related to the escort probability 
distribution through the relation Fq — nPq. The normal- 
ization condition reads J Fqd^U = n. 
We define a collision conserved quantity x(a;, v, t), which 
have the conservation property Xi +X2 = Xfi +X/2- The 
following theorem can be proved: 



fdF\ 

\ dt 
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where is the collisional integral, making the 

changing of variables below: 
First: vi —> V2 and vice- versa. 
Second: vi Vf\, V2 ^ Vf^ and vice- versa. 
Third: vi Vf^, V2 Vf-^ and vice- versa. 
As these changes don't affect the transition matrix el- 
ements Tfi and the delta function, we add each of the 
integrals obtained with equation (lO)and divide the re- 
sult by 4. We end up with: 



d^vx 
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X \TM' {Ff.Ff, - F2F1) (xi +X2- Xh - Xh) 

Through the conservation of xii the theorem is proved. 
Using equation (8), we have, changing the volume ele- 
ment form d^v to d^U: 



d^iXi ( 4 + + — ) Fi(x, v,t)^0 (11) 
ot m 



Defining the expectation value as: 



< O i / FiOd-'U 
n 



(12) 



and rearranging terms in (11), we have the desired con- 
servation theorem: 

d 

— < nxi > +d.t < nviXi > -n < VidiXi > 
ot 



-F, < d^^xi >= 
m 
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Now we calculate the above equation for the three con- 
served quantities, namely mass, momentum and energy. 
For the mass, xi = 'ti, we have: 



where: 



l^-f a,(pu,) = o 



p — mn{x, t) 
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is the mass density of the system. 
For the momentum, xi = fnvi, we get: 

— (pui) + dj {p < UiUj >) + dj {puiUj) = p^ (16) 

Finally for the energy, xi = ^U"^, we have, using the 
mass conservation equation (14): 



d_ n 
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Utilizing parity arguments wc can show that this inte- 
gral will be non-zero only for i = j, so we have, changing 
volume element and the variables as was done in the last 
calculation we obtain: 
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2 /x3/2(l - (jx)T^da; 
m/? / a;i/2(i _o-a;)T^dx 



+ '^diUj < UiUj > +jdi{p < UiU^ >) = (17) 

where we used the fact that the internal energy for this 
system is given by the classical energy equipartition the- 
orem. All we have to do now is calculate the two expecta- 
tion values in order to obtain the hydrodynamic equation 
for the generalized statistics. 



B. The Zero Order Approximation 

To obtain the exact equations, we should use the so- 
lution for (8) Fi. As we don't still have this solution, 
we must approximate it. We assume that the system, 
locally, obeys equation (9). This is the local equilibrium 
assumption. We impose a space and time dependence of 
the macroscopic parameters number density n, average 
velocity Ui and temperature ^. Hence the distribution 
function reads: 
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The integral in the denominator is readily calculated: 

1 



-P—W 



d-'U 



Using the spherical volume element and making the fol- 
lowing change of variable: 
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we have: 
I = 2n\ 



where a is the signal of K-^. It should be noticed that 
this normalization factor depends only on the tempera- 
ture, as in BG statistics. 
We must calculate the quantity < UiUj >q: 



< UiUj > 



g=^j UiUjFg^^U^'U 



For a = +1, these integrals are made in the interval be- 
tween and 1, due to Tsallis cut-off condition, and give 
the result: 



< UiUj >q-- 
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For (T = — 1, the integrals are calculated between and 
-|-oo and give the same result. 

The other expectation value, < UiU'^ >q vanishes in this 

approximation because it is an odd function. 

The mass conservation equation is the same in all orders 

of approximation: 



dp 
dt 
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The momentum conservation equation is written as: 
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Using the mass equation, we have: 
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The term ^ is defined as the pressure of the system: 



PV = NknT 



(22) 



We recognize this equation as the equation of state of the 
ideal gas. As we can see it is independent of q. 
The energy conservation equation is written as: 



d 



P 



-^di{pui)+^5i]^dju. 







dt \(3 

Again using the equation for the mass and dividing by p: 



,1 







(23) 



We conclude that, in the zero order approximation, the 
hydrodynamic equations are q-invariant and we obtained 
the correct equation of state for the ideal gas. 
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C. The First Order Approximation 

To obtain the fully hydrodynamic equations a solution 
of Boltzmann transport equation is needed to be used 
as the distribution function. Since this equation until 
now was not exactly solved, we made an approximation 
in the previous section. Now we add to this zero order 
approximation a correction to it. This procedure consist 
in writing this exact function as: 



- p (0) 



(1) 



where Fq''^^ is the correction to the zero order approxi- 
mation. It can be shown that this correction is given by 
the following relation : 

where r is a parameter of the order of magnitude of the 
collision time. We must calculate each of the derivatives 
of Fq''°\ At first we evaluate the following derivatives. 
For the density we have: 



dp 



d_ 
d~p 



P j(o) — 
ml 



Using these results in the first order approximation func- 
tion we get: 



U,d,pFq^°^ 
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This function is different from the one calculated by the 
Maxwell-Boltzmann function , it has the explicit de- 
pendence on q through the term ^"ttoTi which is unity 

Jig J \ ' 

in the extensive limit. 

The quantity < UtUj >q, in this order of approximation, 
is given by: 



< UiUj >q= - 

n 



(1)^37 



The first term is simply The second term will give 
non- vanishing contribution only through the second term 
from equation (28), the others are odd functions: 
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Second the average velocity: 
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Third the temperature: 
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The derivative of the distribution function with respect 
to the thermal velocity of the molecules is simply the 
derivative with respect to the average velocity with a 
changed signal: 
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This equation is a traceless tensor and it is suffices to 
calculate any component of it. Putting i = I and j = 2, 
we have the result: 

< U1U2 >q= -^rPjdiUk JuiU2UkUi 



1 



Aq ' 2 



This integral is non-zero only for fc = 1,2 and I = 1,2. 
We have, writing Ui = Usen9cos4>, U2 = U sen6sen<j) and 
changing the volume element to a spherical one: 
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Using equations (18) and (19), we have: 
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For the a = +1 case, the above integrals give: 

T 

< U1U2 >q— 3(C^1W2 + d2Ui) 



For the other case, cr = — 1: 



7i = n- 
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The second integral is given by: 
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which results in: 
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(33) 



< C/i?72 >q= ^{diU2 + d2Ui) 

mfj 

So the value of < UiUj >q in this order can be written 
as: 



The last one is given by: 



< UiUj >q= ^ (diUj + djUi - "^SijdkUk 
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This term doesn't change in this context. 

The momentum conservation energy is written as: 



— r (diUj + djUi 

mfj mp 
+ dj{pUiUj) = 
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where the term ^ is the viscosity coefficient of this sys- 
tem: 



mf3 



(31) 



The other expectation value to be calculated is given by, 
where only the first and third terms of (28) will give non- 
vanishing contributions: 



Resulting in: 
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Using the above results, the expectation value < UiU"^ >q 
is given by: 
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The energy conservation equation is given by: 
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Calculating each integral separately, we have: 



U^Fq^°^d^U 



This integral is calculated in the same way that the oth- 
ers above, changing the volume elements, using equation 
(17) and observing the cut-off condition, we obtain: 
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Using the equation for the mass, we have: 
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where: 
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is the q-dependent thermal conductivity coefficient. This 
coefficient, due to its dependence on q, can assume pos- 
itive, g < 1,4 and q > 1,6, or negative values, 1,4 < 
q < 1,6. These negative values are non physical. Also 
we have defined: 



1 — g Srp 
" " A, + (1 - q) 3^ 



T — 



(38) 



which is the new transport coefficient. This one too can 
have positive, q<leg>l,4, as well as negative values, 
1 < (7 < 1,4. This expression vanishes, properly, in the 
extensive limit, g — *■ 1. 

The relation between these two coefficients is given by: 



1-g 

A, 



(39) 



The relation between the thermal conductivity and the 
viscosity coefficients is given by: 



El 
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So the old relation is now q dependent and is equal to 
the old value in the extensive limit since — > 1 when 
q^l. 

The coefficients obtained by Boghosian are related to 
those obtained here through the following relation: 
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Tb 



III. CONCLUSIONS 
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The hydrodynamic equations obtained here are sim- 
ilar to those obtained by Boghosian js). The same 
changes in the energy conservation equation in first or- 
der are observed here, only that we have showed that 
the q-dependence of the thermal conductivity and the 
new transport coefficients are different. We have showed 
also, through the derivation of the conservation equa- 
tions, that the ideal gas state equation is still valid in the 



generalized statistics context. 

As was said there are some differences between our re- 
sults and the ones obtained by Boghosian ||]. The most 
important is that these differences in the g-dependence 
of the transport coefficients, in the thermal conductivity 
and in the new transport coefficient are due the differ- 
ence in the state equation and in the internal energy for 
the ideal gas. In our result we do not have intrinsically 
a g-dependence, while in the Boghosian calculation they 
have a factor Aq = 1 + {1 — q)^. In order to now what 
result is the correct one is necessary to make some appli- 
cation of these transport equations. In despite of this we 
believe that our results are better because we obtain the 
correct state equation (independent of the parameter q) 
and we assume the validity of the classical equipartition 
theorem. 
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